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1. Introduction 

Modern large buildings in the Tokyo metropoli- 
tan area usually have a centralized hot-water sup- 
ply operating 24 h a day. Copper pipe is widely 
used in such systems because of its relatively high 
resistance to corrosion, coupled with additional 
pragmatic merits: it is easy to work with, eas>' to 
install, and relatively cheap. 

The seriousness of Type-II pitting corrosion, 
however, has increasingly received high recognition 
in such hot-water supply systems [1,2,3]. The need 
to obtain information regarding the degree of pit- 
ting corrosion has increased over the last decade 
because considerable pipe damage may require 
maintenance, and even replacement, and in that 
case, proper life prediction is essential to pass rea- 
sonable engineering judgements and thereby to 
perform proper maintenance. 

The life prediction of such copper plumbing tub- 
ing can first be performed by coupling adequate 
nondestructive and/or, though less favorable, de- 
structive inspection techniques with reliable statis- 
tical analysis. It has been shown that the most 



promising statistical analysis methods for such a 
purpose include extreme value statistical analysis 
[3,4]. 

Although such tools have become widespread, a 
general method for evaluating the localized corro- 
sion propensity on existing engineering structures 
from limited inspection data, and concrete crite- 
rion for the number and size of samples required 
to obtain a reasonable extreme value prediction is 
still not available. 

In the present study, a set of pitting corrosion 
depth data obtained from 7 year old copper plumb- 
ing pipe, one third of which was removed from a 
centralized hot-water supply system, was examined 
by extreme value statistics. Emphasis was placed on 
the effect of the total number and size of maximum 
pit depths on the accuracy of pit growth prediction. 
A concept was ultimately proposed to obtain a rea- 
sonable prediction of maximum pit depths while 
minimizing the total sampling area (or, length) for 
analysis. 
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2. Test Procedures 

2.1 Test Specimens 

Copper pipes totaling 8.88 m in length were re- 
moved from various parts of the centralized hot- 
water supply system in an 11-story multi-family 
dwelling In Tokyo. The system consisted of one 
stainless steel storage tank (1,8 m^) and copper 
plumbing pipe having an overall length of about 28 
m. The plumbing material was JIS C 1220T type 
25AM (outside diameter 28.58 mm and wall thick- 
ness 0.89 mm) copper pipe for building use. The 
system had been operating for about 2600 d before 
test piping was removed. In the system, water at a 
nominal temperature of 60 "C circulated con- 
stantly. Average flow (i.e., hot-water consumption) 
was around 8 m^d, the storage tank being supplied 
automatically with tap water. 

The copper piping removed was cut into parts 
100 mm long, which were then cut in half to give 
half-ring specimens. Each half-ring .specimen was 
then completely cleaned ultrasonically in dilute sul- 
furic acid, followed by marking-off to divide it into 
10 virtual half-ring specimens of 10 mm unit length, 
after which the pit depths were measured by using 
an optical microscope of 1 |j,m precision. By cou- 
pling two opposite virtual half-ring specimens, a 10 
mm long full ring specimen was reassembled, and 
the area thus surveyed should be representative of 
the pipe at that particular location. Then, by taking 
several adjoining full-rings of 10 mm length, for 
each unit sample sizes of maximum pit depths 
(hereinafter, s: unit length) in the interval 20 mm 
up to 200 mm were obtained. 

2.2 Extreme Value Statistical Analysis 

The extreme value statistical analysis was per- 
formed by using a commercial available personal 
computer software package, EVANS [5]. 

The basic concept of the present extreme value 
analysis is briefly reviewed in the following sec- 
tions. 

2.2.1 Extreme Value Probability Plots The 
fir^t step of the extreme value analysis included the 
preparation of extreme value probability plots, that 
is, plots of maximum pit depth data on extreme 
value probability paper of the cumulative relative 
frequency (f (y)) vs maximum pit depth (x). Maxi- 
mum pit depths data were arranged in order from 
largest to smallest and assigned a rank number. 
The vertical plotting position F(y) for each pit 
depth value was calculated by the averaged rank 
method ai follows: 



F(y) = }-i/(N + l) 



(1) 



where /= rank number, 
N = total number. 



Several sets of plots were obtained depending on 
the combinations of N and 5. In an exemplifying 
case of s =100 mm and N -7, a total of 10 sets of 
plots was subjected to regression analysis. 

2.2.2 Regression for the Best Fit Line Re- 
gression analysis based on the MVLUE (Minimum 
Variance Linear Unbiased Estimator) method was 
made for each data set to determine a straight line 
of best fit to the plotted extreme values. The equa- 
tion of a straight line is given as follows: 

x = X + oy (2) 

where x = expected maximum pit depth, 
y = standardized variable, 
A =location parameter, 
a=scate parameter. 

The straight line is drawn in Fig. 1 on extreme 
value probability paper. In the present analysis, 
emphasis was placed on the optimum combinations 
of A^ and s for obtaining a reasonable estimate of 
the extreme value. Such an analysis becomes feasi- 
ble through a thorough investigation of all the sam- 
ples where the depth of the actually detected 
deepest pit could be regarded as the probable max- 
imum pit depth, that is, the extreme value. 





, F{y) 




T 




n fti 




100 




U>9J 


1 1 1 


1 1 1 ■ 1 1 


U 










a. 


- 






SO 


at 






/ 




3 






/ 




5 3 


_ 




/ 


20 


3 










E 






/ 




3 
u 


0.90 


- 


/ 


10 


- Z 

o 
u 

c 






A' 


i 


ID 






/ * 


^ 1 


- 0.70 


- 


/ 




u 






f 




u 










o 


0.50 


y 






® 




0.368__ / 






0.30 


f i 






c 










ffi 










3 


0.10 


/ |i 






O" —1 










at ' 










U. 


n ni 


, j . 








Maximum pit depth 

Fig. 1. Extreme value probabiliry paper of double-expunential 
dislribution. 
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3. Results and Discussion 

3.1 Distribution of Pit Depths 

Figure 2 is a histogram showing the relationship 
between the pit depths and the total frequency. In 
the figure, a total of 970 depth data for all pits 
found in one 10 mm long full-ring, removed from 
4he 7th floor in the building, was grouped over the 
pit depth ranges, mm~0.019 mm, 0.020 mm-0.039 
mm, etc. As is evident from the figure, the shape of 
the pit depth distribution is a bell-shaped curve 
starting at zero, rising to a maximum at around 
0.05 mm and thereafter decreasing rapidly with in- 
creasing pit depths. The lack of "J"-shaped portion 
bending to the right in the pit depth range from 
mm to 0.02 mm, together with the trailing extreme 
portion of the tail of the curve up to 0.34 mm, indi- 
cated that most of the small pits had already had 
ceased to grow while only a small number of 
deeper pits continued to grow. 
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Fig. 2. Pit depths distribution histogram fur 10 mm long copper 
pipe used for 2600 d in a hol-water supply system operated ill 
60 °C 

In Fig. 3, all the pit depth data represented in 
Fig. 2 were plotted on a logarithmic-normal distri- 
bution diagram. The apparent linearity of the plot 
indicates that this distribution applies to the loga- 
rithmic-normal; and hence the maximum pit depth 
data obtained in the present pit depth survey were 
sampled from the parent populations with a loga- 
rithmic-normal distribution. 

Figure 4 shows the distribution at every floor of 
the maximum pit depth detected in each 100 mm 
unit length. Though the maximum pit depths 
seemed to have a slight tendency to become shal- 
lower at upper floors, the distribution of pit depths 
was regarded as being uniform throughout the 
building. The actual maximum pit depth value in 
the present survey was 0.452 mm which was de- 
tected at 6th. floor. 
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Fig. 3. Plot on a logarithmic-normal distribution diagram of all 
data presented in Fig- 2. 

In Figure 5, the maximum pit depth data for unit 
lenghs of 20 mm, 100 mm, and 200 mm were evalu- 
ated from extreme value analysis on the basis of 
Gumbel's double-exponential distribution; the lin- 
earity of each plot shows that this distribution ap- 
plies to the maximum pit depth data obtained at 
the unit lengths between 20 mm and 200 mm. 

3.2 Minimum Required A'^ and s 

In the practical application of extreme value 
statistics, the number and the size of unit samples 
for the pit depths survey are to be decided prior to 
the destructive (or nondestructive) inspection of 
the existing structures. 
Reliability of the extreme value prediction de- 
pends on the following three equations; 



>'=lnr 
T=S/s 



(4) 



V(x) = a'\A(N^)y'+B(N/i)y + C(N^)] (5) 



323 



Volume 99, Number 4, July-Augusi 1994 

Journal of Research of the National Institute of Standards and Technology 



.a 

E 



o 

o 






_L 



_L 



Q 0.1 0.2 0.3 0.4 0.5 
Maximum pit depth (mm) 

Fig. 4. Distribution at each floor in the building of the maxi- 
mum pit depth detected in each 100 mm unit length. 



where y = standardized variable, 
7*= return period, 
S - total length, 
i =unit length, 
V{x) = variance of the estimated extreme 

value, 
A(N,n), B(N,n), C(/V,«) = MVLUE 

coefficients, 
a = scale parameter. 

It is obvious from these equations thai for in- 
creased reliability, hence for a minimized K(jc), it is 
required to increase s (therefore to decrease 7") 
and/or to increase N. In the interest of economy, 
however, there must be a natural limitation to the 
increase of N and s for the extreme value survey. 

In theory, the optimum combinations of N and s 
can be determined from Eqs. (3) to (5) once the 
distribution parameter (that is, a ratio of the loca- 
tion parameter to the scale parameter, o/A) is 
reasonably assumed (that is, empirically or experi- 
mentally), and the extent of the standard deviation 
of the error of extreme value estimates, a, may be 
expected to be: 



A =mcr 



(6) 



where A = location parameter, 

(r= standard deviation of the estimated 

extreme value, 
m = assumed number (1, 2, 3, etc.) 

That is, the decision may be made to reduce a to 
1/m of A [6]. 
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Fig. S. Double-exponential distribution t)pe of extreme value probability plots for j = 20 mm, 100 mm, and 200 mm. 
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Table 1 shows the result of extreme value analy- 
sis to obtain the relationship between unit length 
and distribution parameter. In the analysis, all pit 
depth data obtained from the whole lengths were 
brought into consideration. It can be seen in the 
table that the distribution parameter, a/A, de- 
creased with increasing unit length, ultimately ap- 
proaching to a definite level of 0.15. 

Tabk 1. Location, scale, and distribulion parameters deter- 
mined from alt maximum pit depth data obtained in the survey 



Unit length 
(mm) 


T 


N 


(mm) 


A 

(mm) 


a 

(mm) 


a/A 


50 
100 
200 


177.6 
88.8 
44.4 


172 
36 

43 


0.575 
0J28 
0J06 


0.28415 
0J0446 
0J221 


0.05621 
0.04986 
0.04866 


0.198 

0.16 

0.151 



Optimum combinations of N and T (hence, s) 
required for controlling o- down to 1/m of A were 
obtained as shown in Fig. 6 corresponding to a/A of 
0.15 and 0.20. It stands to reason, that the number 
of unit lengths, A^, can markedly be reduced by 
loosening the requirement for the reliability, that 
is, by decreasing m . N may also be reduced by de- 
creasing a/A. Though N may be reduced by in- 
creasing s (hence decreasing T), it should be noted 
that the total length required for a pit depths sur- 
vey can increase on the contrary. 




Ftg. 6. Optimum combinations of l^ and T to control a at levels 
of 1/2, 1/3, 1/4. or 1/5 of A under the limitation of tt/A =0.2. 



3.3 Maximum Pit Depths Prediction Based on the 
Optimized N and s 

Figure 7 shows the results of the maximum pit 
depths prediction by using the combinations of N 
and s from Fig. 6. Unit lengths for this evaluation 
were not sampled at random locations, but were 
ordered from one end of the lower floors upward. 
Thus the results were represented by correlating to 
the locations wherefrom those unit lengths used in 
the prediction were removed. It can be seen that 
the scatter of estimated extreme values decreased 
with increasing s. As would be expected, the scat- 
ter was narrower at m =5 as compared with m =3. 

To determine the general tendency in the reli- 
ability of the present extreme value prediction, 
Table 2 was developed from data on Fig. 7 by tak- 
ing an average for each item. The results indicate 
the following: 

1) The assumptions for a/A and m, that were 
made prior to the analysis, have conservatively 
been met. 

2) Increased m did not always result in increased 
reliability which indicated that m =2 might be 
reasonable in the interest of economy. 

3) Except for the cases of 5 = 10 mm, the maximum 
detected pit depth value of 0.452 mm fell closely 
between the average of estimates ±la-. 

4) Although the level of s is to be kept as low as 
possible, because the total sampling length (ie., 
N xs) may increase with increasing s, at least 
2.5% of the whole structure is to be subjected to 
extreme value evaluation at a return period of 
500 or the less. 



4. Concluding Remarks 

Based on this analysis, the following conclusions 
can be drawn concerning the optimum conditions 
for obtaining a reasonable extreme value predic- 
tion. 

1) The number and the size of unit samples for the 
extreme values survey may be determined so 
that the variance of the extreme value estimates 
is to be minimized under a definite distribution 
parameter, wherein the standard deviation of 
the estimates is expected to be 1/3 of the mode 
of distribution. 

2) A return period of 500 with a total sampling 
length amounting to 2.5% of the entire parts 
may be desirable for an increased reliability of 
extreme value prediction. 
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Fig, 7, Distribution al each floor of the estimates of maximum pit depths on the hasis of the optimum 
combinations of W anil 7" under the limituliunx of A =3o-and 5v. 

Table 2, Summary of the extreme value prediction 















Estimated 




A„ 




Unit length 


Return 


Sample 


Location 


Scale 


maximum 


OXaat 




S 


period 


number 


parameter 


parameter 


pit depth 




•ifitm 




(msl) 


r 


A' 


A„ 


Qjv 


Jfm. 








10 


888 


14 


0.241 


0.043 


0.534 


0.069 


3.50 


A=3ff 


20 
40 


444 

222 


12 

9 


0.269 
0.288 


0.034 
0.035 


0.478 
0.467 


0.054 
0.058 


4.98 
4.98 




100 


88.8 


7 


0J13 


0.028 


0.428 


0.045 


6.94 




10 


888 


33 


0.218 


0,046 


0.559 


0.047 


4.63 


A=5o- 


20 


444 


2B 


0.266 


0.040 


0.491 


0.040 


6.67 


40 


222 


23 


0.282 


0.039 


0.487 


0.039 


7,25 




100 


8S.8 


18 


0.312 


0.034 


0.463 


0.032 


9.71 
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